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We have used existing data on the one-proton-removal cross section of at 285 MeV/u and 
Glauber model calculations to extract the asymptotic normalization coefficient for the wave function 
of the last proton in the ground state of ^C. The calculations are done first using folded potentials 
starting from two different effective nucleon-nucleon interactions and second in the optical limit 
using three nucleon-nucleon interactions, and the results are found to be consistent, with no new 
parameters adjusted. We find C^{p^/2) + C^{pi/2) = 1-22 ± 0.13 fm~^. From this result we obtain 
the astrophysical factor for the proton radiative capture reaction *B(p,7)^C as 518(0) — 46±6 eV-b. 
' The calculated energy dependence of the astrophysical S-factor for the energy region Ecm = — 0.8 

MeV and the reaction rates for Tg = — 1 are included. 
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I. INTRODUCTION 

(N : 

>• . . . 

, Radiative capture reactions are of crucial importance in nuclear astrophysics. The capture of charged particles 
on nuclei, in particular of protons in basic processes like hydrogen burning in various stellar environments, is very 
much hindered by the Coulomb repulsion. This leads to very small cross sections and consequently to the well 
known experimental problems associated with direct determination of astrophysical S'-factors at very low energies 
0. Moreover, with the large number of reaction chains found to be of importance in nucleosynthesis calculations for 
different static and explosive burning scenarios l^j^Q, more data involving the capture on unstable nuclei becomes 

' necessary. In many instances direct measurements involving unstable nuclei are very difficult or even impossible and 
^ indirect methods must be used. 

1^ It has been known for a long time that transfer reactions can be used as indirect methods for nuclear astrophysics 

(see, e.g., and references therein) and various techniques have been used since. A few years ago another indirect 
approach based on measurements of peripheral transfer reactions was proposed , and was subsequently used to 
determine the astrophysical factor 6*17(0) for the reaction ^Be(p,7)^B using (^Be,7B) transfer reactions. The method, 
which is based on the observation that radiative capture of protons is a very peripheral process at astrophysical 
energies, involves the extraction of nuclear quantities called Asymptotic Normalization Coefficients (ANC) from 
proton transfer reactions and these ANCs are then used to determine the astrophysical cross sections. The method 
5J] works if the transfer reactions, which have much higher cross sections than the capture reactions, are also peripheral. 
Ci ^ Using secondary beams available today, such experiments can be done in a matter of days [pp|, p^ , pl] , p^ . More 
recently we have shown fl^ that one-nucleon-removal reactions offer an alternative and complementary technique for 
extracting ANCs which is particularly well adapted to rare isotopes produced using fragmentation. In the breakup 
of loosely bound nuclei at intermediate energies, the requirement of core survival in the final channel automatically 
selects core-target interactions which are highly peripheral, with the result that the wave function of the removed 
nucleon is probed at and beyond the core's nuclear surface. This is then very similar to the low-energy light ion 
transfer reactions where the short mean free path of the ions in the optical potential leads to surface localization. 
In Ref. *B breakup data at energies between 30 and 300 MeV/u on a variety of targets were compared with 
calculations done using a Glauber model to extract the same ANC as that obtained from transfer reactions at 
about 12 MeV/u. We suggested that such experiments have the additional advantage of using beams of lesser quality, 
at higher energies, and therefore can be applied to nuclei farther away from the stability. 

Here we use a similar approach to analyze the breakup of ®C. Existing cross section data for the one-proton-removal 
reaction with ^C at 285 MeV/u on four targets (l5| are used to extract the ANC for the virtual process ^C— >*B-)-p. 
The breakup cross sections are calculated with a Glauber type model, first using the potential approach described 
in now with two types of effective interactions, and then using a Glauber model in the optical limit, with three 
different prescriptions for the elementary nucleon-nucleon (NN) scattering amplitudes. Good consistency is obtained 
between the different approaches. The value of the ANC is then used to determine the astrophysical factor Sis for 
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the proton radiative capture reaction ®B(p, 7)^0 in the energy range Ecm = — 0.8 MeV and the reaction rate in 
stellar environments is obtained for the temperature range Tg = — 1. The ANC extracted here agrees with that 
obtained recently from a transfer reaction d(^B,^C)n on a deuteron target at 14.4 MeV/u but our result has a 
smaller uncertainty. 



II. THE BREAKUP OF '^C 



Much effort has been devoted in recent years to the study of exotic nuclei. In particular, it has been shown that 
onc-miclcon removal reactions with radioactive projectiles at intermediate energies can be used to study their structure 



[ [16| , |17 18 19[]. Typically an exotic nucleus B = (Ap), where B is a bound state of the core A and the nucleon p, is 
produced by fragmentation from a primary beam, separated and then used to bombard a secondary target. After 
breakup occurs, the cross section and parallel momentum distribution of the core A are measured. It has been shown 
that the parallel momentum distribution of the detected core A can be used to establish the orbital momentum 
of the relative motion of the nucleon, even disentangle contributions from different contributing orbitals, and that 
coincidences with gamma-rays enables one to separate the contribution of different core states into the ground state 
of the projectile. Typically the integrated cross sections have been used to extract absolute spectroscopic factors. 
However, we have shown [ p"3|j20| ] that the extracted spectroscopic factors depend strongly on the assumed single 
particle wave function (or cquivalcntly on the geometry of the binding two-body potential that produces it), and 
unfortunately, for exotic nuclei these potentials are poorly known and there is not much other information to put 
constraints on them. We have shown that because these breakup processes are essentially peripheral, especially for 
loosely bound nuclei, one can obtain the asymptotic normalization coefficients (ANCs), rather than the spectroscopic 
factors. The ANC only gives information about the asymptotic tail of the wave function for the last nucleon, but does 
not say anything about the behavior of the many body wave function inside the nucleus. Fortunately, this is all we 
need to determine peripheral observables, chief among them the astrophysical S'-factor for radiative proton capture 
reactions. The extraction of ANC also has the advantage of not depending on the fundamental assumption about the 
validity of the independent particle model in the whole nucleus, but only at its surface (see |^,^ and the references 
therein, for a discussion on this subject), or of the cluster-like structure of some loosely bound nuclei, an important 
characteristic. 

We apply here the method described first in |p^ , to analyze existing one-proton-removal cross section data for ^C 
at 285 MeV/u on four different targets (C, Al, Sn and Pb) ||l^. In these experiments a beam of radioactive ^C nuclei 
stroke the target, and the residual *B was detected. We use the same extended Glauber model ^3 24,2^] to calculate 



the cross sections. The applicability of such calculations for breakup reactions was discussed before in Refs. | |26| , p 
In all reactions where the core survives (either proton transfer or one-proton breakup) the matrix elements for the 
transition B ^ A + p include the overlap integral I^p{r) for the nuclei A, p, and B, obtained after integration over 
the internal coordinates of A for the many-body fully antisymmetric wave functions, with r the vector connecting the 
center of mass of nucleus A with p . The ground states of the loosely bound nuclei are known to be dominated 
by single particle features. Outside the core both the nucleon-nucleon correlation effects and the antisymmetrization 
effects are small and the overlap integral behaves very much like the radial wave function for a single particle in the 
potential given by the core : 

jB (r)^ s'/^w ,-(r)^r^ W^-^,i+i/2(2^r) 

Here Snij is the spectroscopic factor and, in the rightmost part of the equation, C^^ij is the asymptotic normalization 
coefficient defining the amplitude of the tail of the overlap integral, W is the Whittaker function, n is the wave number, 
and 77 is the Sommerfeld parameter for the bound state (Ap). Due to Coulomb repulsion, the radiative proton capture 
at stellar energies takes place at very large distances from the core, in regions where approximation (Q) is very good [Q, 
but the overlap integral there, and consequently the resulting cross sections, are very small. However, the asymptotic 
normalization coefficients C^pij be extracted from phenomena that are peripheral, but involve regions of radii 
much closer to the core, and therefore have larger cross sections and are easier to measure experimentally. Another 
important reason for larger cross sections can be, of course, the different nature of the transition operators involved. 
One-proton removal reaction from loosely bound projectiles is one such phenomenon. As an example to illustrate the 
above statement, for the well known reaction ^Be(p,7)^B at solar energies (Ecm=20 keV) the classical turning point 
in the entrance channel is around r « 250 fm, but the combined effect of barrier penetration, transition operator and 
asymptotic behavior of the ^B bound state wave function makes that regions around r « 50 fm contribute most in 
the radiative capture process. This means that the wave function of *B is most effectively sampled at these distances, 
where it is very small. Going to processes where distances r w 4 — 5 fm contribute most, like in the case of breakup 
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(fig.l in Ref. |l^]), one can gain 5-6 orders of magnitude in cross section from the magnitude of the wave function 
alone, as suggested quahtatively above, while the only unknown quantity is the same ANC. We treat the breakup 
reactions in this section using two approaches in a Glauber model. 

A. Glauber model calculations with folded potentials 

In the extended Glauber model used here, the center of mass of the projectile (made up of core and loosely bound 
proton) moves on a straight line trajectory, an approximation valid at intermediate energies. The relative motion of 
the outer proton about the core is described by a single particle wave function with the asymptotic behavior given 
by Eq. (|^). They are assumed to interact separately with the target. The breakup of the projectile appears from 
three different processes: the proton is absorbed by the target while the remaining core is scattered and detected 
(stripping) ; both the proton and the core are scattered by nuclear interaction with the target and the core is detected 
(diffraction dissociation); dissociation in the Coulomb field of the target. The probability for each process depends on 
the impact parameter of each trajectory and on the wave function for the proton-core relative motion. The total cross 
section for this one-proton-removal reaction is given by the integral over all impact parameters. S-matrix elements 
have been calculated in the eikonal approximation including corrections up to second order to assure convergence. 
The S-matrix is given in this potential approximation by 

S{b) = e'^^''\ (2) 

where the leading term in the eikonal expansion is calculated along the trajectory 

X{b)^-^ r dzV{b,z). (3) 

The S-matrix elements are calculated for each trajectory, the contributions to the three terms above are evaluated 
and then summed incoherently as 

(7sp ^ (Tstr + <7diff + CrCoul- (4) 

For the reaction model calculations we assume that the ground state of the projectile {J'") can be approximated 
by a superposition of configurations of the form [I^" (g) nljY , where IJ!" denote the core states and nlj are the 
quantum numbers for the single particle wave function Lpnij{r) in a spherical mean field potential. These single 
particle states are normalized to unity and have the asymptotic behavior given by Eq. with the single particle 
asymptotic normalization coefficients When more than one configuration contributes to a given core state, 

then the total cross section for one-nucleon breakup is written as an incoherent superposition of single particle cross 
sections, weighted by the spectroscopic factors S{c,nlj): 

o-_ip = ^ S{c, nlj)asp{nlj). (5) 

For ^C the last proton mainly occupies the lp3/2 and lpi/2 orbitals around a core (which only has one bound 
state), in unknown proportions. At large distances the two orbitals have identical radial behavior and consequently 
have the same single particle cross section asp and the same single particle asymptotic normalization coefficient bp. 
The one-proton-removal cross section (^ can be written as the sum of two components: 

= [S{lpy,) + S{lp,/,)] asp{lp^) = {Cl^^ + Cl,,)osp{lp,)lhl (6) 

and we can only extract the sum of the two spectroscopic factors or the sum of the two (squares of) ANCs C^ff — 
^P3/2 ~^ This same sum appears in the evaluation of the astrophysical S'-factor. Here we have only considered 

the lp3/2 orbital. This was taken as the solution of a radial Schrodinger equation in a Woods-Saxon potential with 
a geometry given by its half-radius and diffuseness parameters {R,a) and a spin-orbit term with strength Vl5=18.6 
MeV-fm^ pi. The depth of the potential was adjusted to reproduce the experimental separation energy of the last 
proton Sp = 1.296 MeV, for each geometry (different values of {R, a) used). 

The coordinate system used in the Glauber calculations is shown in Fig. 1. To calculate the S-matrix elements we 
need the interaction potentials between the three participants. For the proton-target potential we have used, as we 
did before, the G-matrix effective interaction of Jeukenne, Lejeune and Mahaux (JLM) in the updated version 
of Ref. Q. This interaction is complex, energy and density dependent and was adjusted for a very wide range of 
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targets and proton energies. For the target-core potential the double folding procedure described in was used: the 
same JLM interaction was folded with Hartree-Fock nuclear matter distributions of the core and of the target. The 
single particle densities were obtained in a spherical Hartree-Fock calculation using the density functional method 
of Beiner and Lombard The strength of the surface term in the functional was adjusted slightly in order to 

reproduce the known experimental binding energy for each nucleus. This procedure, similar to those used in Refs. 
1^,^, leads to root-mean-square radii of the proton and mass distributions < >^l^—1.hl fm and < r^^ >^/^=2.43 
fm, respectively, for ^B, in fair agreement with our experimental determination using the ANC method ||3^ . The local 
density approximation was improved to include finite range effects by using smearing normalized Gaussian functions 
of ranges Ib. and t/. The resulting double folding potentials were subsequently renormalized to reproduce a variety 
of elastic scattering data for light nuclei. We found in Ref. that at incident energies of about 10 MeV/u the 
real part of the potential needed a smearing with a Gaussian of range = 1.2 fm and a substantial renormalization 
{Nv — 0.37), while the imaginary part did not need renormalization {N\y = 1.00) and the smearing was larger 
tj — 1.75 fm. In the present calculations we adopted this procedure, with the JLM(l) interaction, and Nw — 1.00. 
The procedure and its parameters determined by this approach were found to give good predictions for the elastic 
scattering of radioactive nuclei such as ^Be [|j and ^^C ||l^ at similar energies. We have also checked the procedure 
on a much wider set of data from the literature at larger energies, and found a similar conclusion for the imaginary 
part of the potential, while the renormalization of the real part approaches unity around 50 MeV/u. The S-matrix 
calculations that enter the first two terms of Eq. ^ depend primarily on the absorption, and thus on the geometry 
and strength of the imaginary part of the interaction. The real part of the potential only influences the phases of 
the scattered waves as the particles go along the classical trajectories and does not affect the total cross section. The 
Coulomb dissociation term is treated in a perturbative method equivalent to that of Ref. | |39[ |, except that radial 
matrix elements for El and E2 transitions were calculated with realistic Woods-Saxon radial wave functions. Thus 
both nuclear and Coulomb breakup contributions were calculated in a consistent way. 

In order to check that the process was peripheral, the stripping and the diffraction dissociation probabilities were 
calculated as a function of the parameter s (Fig. 1), which is the projection of the proton position vector relative to the 
core, onto the plane perpendicular to the beam direction. The stripping and the diffraction dissociation probability 
distributions are shown in Fig. 2. It is clear that for all four targets (the calculations were done for single isotopic 
targets ^^C, ^^Al, ^^^Sn and ^°*Pb) the contributions from both terms peak outside the radius of the assumed ®B 
core, Rc=2.60 fm However, the interior does contribute and it should not be discarded. Note that for all four 
targets the stripping cross section dominates, with the diffraction dissociation part contributing only about 15%. This 
situation is different from that depicted in Fig. 1 of Ref. ||l^] where the same two probabilities are shown for the 
case of ®B breakup on a Si target at 38 MeV/u. In that case the two components of the total cross section have 
similar magnitudes, with the diffraction dissociation part being slightly larger. This is a well understood bombarding 
energy effect. High energy heavy ion interactions are dominated by strong absorption. The ®B calculations are in 
remarkable agreement with the experimental data that explicitly disentangled the stripping and dissociation reaction 
mechanisms 123]. The present results are for a much larger energy, and the stripping term (absorption) dominates, 
as it does for°B breakup at larger energies. The third term that contributes to the total one-proton-removal cross 
section. Coulomb dissociation (acoui), is the most peripheral due to the long range of the Coulomb forces. It was 
calculated in its integral form and is not included in Fig. 2. In Fig. 9 of Ref. |4^] its radial dependence is shown 
for situations close to the present one. The present calculations were repeated, changing only the half-radius and the 
diffuseness parameters of the proton binding Woods-Saxon potential that describes the relative motion between the 
proton and the core on a 20 point grid using R = 2.20 — 2.60 fm in 0.10 fm steps and a = 0.50, 0.60, 0.65 and 0.70 
fm, but keeping the reaction mechanism parameters unchanged. 

From the comparison of the calculations with the experimental cross sections we extracted C^ff — , using 

Eq. g. The results are compared in Fig. 3, in order, left to right, for the C, Al, Sn and Pb targets. The error bars 
contain the contribution of the experimental errors and of our uncertainties in the calculations, which are described 
in detail below, added quadratically. The distribution of the four numbers extracted is consistent with the constant 
value C^ff — 1.18 ± 0.12 fm~^. The results are also shown in the third column of Table I. 

In a second approach, we used the free NN i-matrix interaction of Franey and Love |^l| in a tpp approximation to 
obtain the interaction potentials. This is a local representation of the free NN interaction based on phenomenological 
nucleon-nucleon scattering amplitudes of Arndt et al. at several energies between 50 and 1000 MeV. In each 
spin-isospin channel the interaction is given by a linear combination of Yukawa form factors. The longest range is 
fixed to be the long range part of the one pion exchange potential (OPEP), but this part does not survive in the 
direct term. In the calculation of nucleon-target and projectile-target interactions, the free NN t-matrix should be 
renormalized. We used the prescription suggested in pl[ employing relativistic kinematics. At intermediate energies 
one should interpolate in the Franey-Love tabulation. In practice, this procedure is difficult since the number of 
components varies from energy to energy. At 285 MeV/u we used the parametrization at the nearest energy (300 
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MeV). Other effects such, as Fermi averaging, Pauh blocking and spin-orbit contributions were ignored. Finally, the 
scattering phase was calculated with Eq. (g) in the lowest order of the eikonal approximation. We have checked 
however that higher order noneikonal corrections are negligibly small at this energy. The values extracted for C^ff 
are similar with those obtained with the JLM interaction and are presented in the forth column of Table I. 



B. Glauber model calculations in the optical limit 

At 285 MeV/u, the optical limit of the Glauber model is also applicable. In order to check the cross sections 
calculated above in the folding model, the S-matrix elements were generated in the optical limit which is just the 
leading term of the cumulant expansion of Glauber's multiple scattering theory: 

Xib) = ^crNN{i + unn) J dbidb2Pprojibi)ptarg{b2)v{b + bi- h), (7) 

where Pproj{targ){b) is the profile density of the projectile (target) obtained from the same Hartree-Fock densities used 
before: 

/ + 00 
dzpHF{b,z) (8) 
-oo 

(z is the beam direction) and a^N is the ratio of the real to the imaginary nucleon-nucleon forward amplitude. These 
were interpolated for each isospin state from Table I of Ray gnn is the NN total cross section at the given 

energy, taken from the parametrization of Ref . . The integral appearing in Eq. is a projection onto the impact 
parameter plane of the interaction potential 

ri{R) = j dfidf2Pproj{ri)ptarg{r2)v{ri+ R- r2) (9) 

Following Ray we identify the Fourier transform of the interaction v with the elementary scattering amplitude: 

Hq) = —f = YCTjvw(awjv + z)e~''""« (10) 

where q is the transferred momentum and v is the asymptotic velocity. In practical terms Eq. (|^) was solved using a 
suitably chosen elementary interaction v, such as 

dre*'~w(r) = e"^'^^ (11) 



with the solution 

This was projected onto the impact parameter space and using Eq. (^ the scattering phase was calculated. The 
amplitudes unn and the ranges p = \/4fi were interpolated from Table I of Ray ji^. In the zero range approximation, 
/i — > 0, Eq. (^ is just the overlap volume of the two densities and the scattering phase is determined by the elementary 
reactions in this volume. The breakup cross section is highly sensitive to the range parameter. To assess this sensitivity 
we have performed calculations with the ranges prescribed by Ray, with zero-range for all interactions and then with 
an average value of /i = 1.5 fm in all pp, nn and pn channels. This last approach is denoted below as the "standard" 
procedure, as it is widely used to estimate the size of halo nuclei from reaction cross sections [Q. When the incident 
projectile was a nucleon, then the corresponding profile density was taken as a S function. We use the charge 
independence assumption i.e. app = ann, cTpn = cFnp, ctpp = ctnn, otpn = anp- 

One test of the model and of the parameters used to calculate the breakup cross sections, was to see how they 
reproduce the cross section data for the assumed components of the projectile {p and ®B), separately, on a ^^C 
target (where the nuclear component dominates). We have calculated the reaction cross section and the total cross 
section for p+^'^C for energies E=100 to 1000 MeV, and the reaction cross section for and for in the energy 
range E/A=100 to 1000 MeV/u. The results are compared with experimental data in Fig. 4. The proton reaction 
cross section data (open circles) in the top panel are taken from the NNDC data base [Ea and those for the total 
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cross section (full circles) are from Ref. [Q. The reaction cross section data for ®B shown in the middle panel are 
from Ref. (filled circles) and Ref. (open circle), while the reaction cross sections for in the lower panel are 
again from Ref. (filled circle) and from Ref. (open circles) , respectively. The calculations shown were done in 
the zero-range approach, and give a good description of the data over the whole range of energies. Therefore we have 
confidence that we can proceed to breakup cross section calculations. We did the calculations with the three range 
parametrizations described above: Ray, standard and zero-range. For each paramctrization, the calculations were 
repeated on the grid of (R, a) parameters for the proton binding potential, as described above for the JLM potential 
calculations. Further details of all these calculations will be presented elsewhere ]25| ]. The results are summarized 
in Table I and are shown in Fig. 5, where they are compared directly with the values already shown in Fig. 3. 
The error bars are the same as in the preceding figure, and the different points show the values obtained using the 
three different prescriptions, as explained above. A good consistence can be observed. We underline that the two 
Glauber approaches are different in essence, and no further parameters were adjusted here. The calculations done in 
the zero-range approximation gave smaller cross sections (and consequently larger ANCs result) by 5% to 16% (less 
for Pb, more for ^^C target) and we exclude them from the overall average, as the zero-range prescription is not very 
realistic. A similar check, done with a larger range of the NN forces (/i = 2.5 fm) shows that the extracted ANCs 
decrease by 10 to 25%. These also are not included in the average. These two extremes in the range parameter provide 
a lower and an upper limit on the cross sections. The final average includes the results obtained with four sets of 
calculations: JLM, Franey-Love, "standard" and "Ray" NN interaction. The weighted average is C^jf = 1.22 ± 0.13 
fm^^, which agrees well with that found from the calculations in the eikonal approach using the JLM interaction. The 
uncertainty takes into account the experimental uncertainties, and the uncertainties in the calculations. The ANC 
values were first averaged for each target, and then an average was calculated for the four targets, using weights given 
by the quadratically combined experimental and calculational uncertainties. To evaluate the uncertainties for the 
calculated values we considered independently those arising from the assumed nuclear structure (the geometry of the 
proton binding potential), as described below, and those from the reaction model. For each target we assumed that 
the standard deviation of the values obtained from the four different calculations give a fair measure of the reaction 
model contribution to uncertainty. Averaging first for each type of calculations (as in the last row in Table I), and 
then for all cases, leads to similar numbers {C^jf = 1.24 ±0.13 fm~^), in part because the experimental uncertainties 
dominate. 

A discussion of the validity of the approach used here and of the uncertainties occurring for the particular case of 
at this energy, a posteriori, is in order. The use of such a fragile nucleus as is as a core can be questioned, and 
also the high energy (285 MeV/u) may lead to questions about the peripheral nature of the reaction. These questions 
are supported by the non-negligible contributions of the interior of the nucleus, as shown in Fig. 2. The proton 
separation energy is not as small here (1.296 MeV) as it was for '^B (0.137 MeV). Using the procedure of Ref. [^ , 
and the ANC extracted we find the rms radius of the last proton in ^C to be (r^)^/^ = 3.02 fm, considerably smaller 
than that found for ^B: (r^)^/^ — 3.97 fm. This, and the large separation energy, reflects in a less sharp definition 
of the extracted ANC. In a figure equivalent to Fig. 2 of Ref. ||l^, the ANC line has a larger slope, pointing toward 
the conclusion that there is a residual dependence on the geometry of the proton binding Woods-Saxon potential 
chosen. This contribution to the uncertainty varies between 7.8% and 4.6% for the targets between C and Pb, and 
is clearly smaller than the experimental uncertainties, and was considered in the uncertainties plotted in Figs. 3 and 
5. No new parameters were adjusted in the present calculations. However some parameters were extracted before, 
with corresponding uncertainties. For the case of the extended Glauber model calculations with the JLM interaction, 
these parameters are the renormalization (A^_r, N\y) and the finite range {Ir, ti) parameters of the real and imaginary 
parts of the nucleus-nucleus potential and have uncertainties. To see how these propagate to calculated breakup cross 
sections, we repeated the calculations varying individually each parameter. We discussed above that the real potential 
does not infiuence the total cross section. Thus only the influence of ti and Nw was determined. We found the relative 
variation of the extracted ANC ^Cg^j/Cg^j to be from ±6.4% on the C and Al targets, down to ±3.3% on the Pb 
target (the relative contribution of the nuclear breakup decreases as the charge of the target increases), when the 
renormalization coefficient Nw varies in the range found in Ref. Nw = 1-00 ± 0.09. Similarly, a change of 4.7% 
to 3.0% was found for a variation of 10% around the central value of the smearing range ti = 1.75 fm. Combining the 
two in quadrature we obtain an 8% uncertainty, which is smaller than the experimental uncertainties are, and similar 
to the uncertainties from the choice of the model. For the calculations in the optical limit, the differences between the 
three approaches used give an assessment of their overall uncertainty. The maximum deviation around the average is 
less than 7-8% in each case. Our present model does not consider two-step processes like ^C-^^B+p ^Be +p + p 
and we cannot therefore directly assess their importance. However, a recent paper ||49[| compares Glauber model 
calculations using a three-body projectile model with those using a two-body projectile (as used in our calculations) 
and finds the dynamical correlations to have a small influence (less than 10%) on the calculated one-neutron-removal 
cross section for the (^^Be,^^Be) and ('^^C,^^C) reactions, a finding that supports the spectator assumption even when 
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the spectator itself is a loosely bound system, as in the case of ^C. The situations considered there are similar to the 
present one, as ^^Be (assumed core for ^^Be) and ^^C (core for ^^C) are halo nuclei of an even larger extent than ^B 
(core for ^C) and, therefore, expect the effect of two-step processes to be even smaller for our case. 

In a recent publication jl^, the ANC for ^C^^B+p was found using the proton transfer reaction d{^B,^C)n at 14.4 
MeV/u incident energy. The experimental statistics were rather poor, according to their Figure 1, and the authors 
present the results of a range of DWBA calculations, with different optical potentials from literature. They report an 
effective ANC that ranges from 0.97 to 1.42 fm~^, with an average that we found to be {C^ff) = 1.18 fm^^. From 
their assessment of the final uncertainty we figure ^C^ff = 0.34 fm~^. This uncertainty is large (30%), and {d,n) 
reactions have been criticized before |Q for not being good peripheral reactions, inadequate for the determination of 
the ANC. However, their values (and particularly the average) are very close to those extracted by us from different 
experimental data and using a different reaction mechanism. 



III. THE ASTROPHYSICAL FACTOR Sis 

Since the ^B(p, 7)^C capture process at astrophysical energies is a peripheral process, the absolute normalization of 
its astrophysical 5*18 factor is entirely defined by the ANC for ^B +p ^C. To calculate the astrophysical S'-factor we 
used the potential model, as described in Ref. |5^. Electric dipole and quadrupole transitions were included for the 
final channel, with El giving the largest contribution, and practically all waves were considered in the entrance channel 
(but the s-wave dominates the major El term and the d-wave contributes only a few percent). The calculations were 
done with a single proton lp3/2 wave function normalized to unity and having the asymptotic normalization coefficient 
bp. Then the result was scaled by Cl^f/h"^ (such a procedure avoids any complications that might appear when a 
Whittaker function normalized by Ce// is used over the whole integration range). The calculations were done for the 
proton energy range Ecm = — 0.8 MeV. The contribution of the resonant state at Eres — 922 keV with known width 
r = 100 keV was found to be unimportant here, because is rather far away and most likely its spin is = 1/2^ 
and thus it is forbidden by selection rules to contribute to the major term. A very weak dependence on energy is 
observed: S{E) = 45.8 - 15.1 • E + 7.34 • E'^ {E in MeV), less than 15% decrease on the whole range. 

Therefore, we find a value 5'i8(0) = 46 ± 6 eV-b. This value is in agreement with the value 518(0) = 45 ± 13 eV-b 
reported in Ref. ||ll[, but disagrees with those previously obtained from various calculations. This radiative capture 
process was investigated in a microscopic cluster model using the Generator Coordinate Method with two different 
potentials (Volkov V2 and Minnesota) in Ref. |5^. The author obtained for 5i8(0) values around 80 and 90 eV-b, 
respectively, therefore 80 — 100% larger than our value. Earlier, Ref. Q used a single particle model to estimate 
the S-factor and obtained a much larger value (210 eV-b), almost certainly because of the too large spectroscopic 
factor C^S = 2.5 that they use. A short communication by Timofeyuk (53) reported a value 518(0) = 53 eV-b 
from microscopic calculations using the M3Y nucleon-nucleon interaction, via an estimate of the ANC and of the 
equivalent vertex constant, but without giving too many details about the calculations. This last calculation is closer 
to the experimental value determined here. Significant disagreement between the astrophysical factors obtained in 
p2| and to a certain extent in |Q using different potentials (although in slightly different microscopic approaches) 
demonstrates the sensitivity of the calculated ANCs to the effective NN-potentials adopted for the nuclear structure 
calculations. Such sensitivity has been observed before for the Ip-sheW nuclei in Ref. |Q. Later the Generator 
Coordinate Method (GCM) calculations with Volkov potentials led to an ANC for the system ^Be +p — > ^B 
about 70% higher than the ANC determined experimentally from proton transfer reactions |^ and from the breakup 
reaction of *B The most instructive dependence of the microscopically calculated ANCs on the adopted effective 
NN-potential has been demonstrated by Baye and Timofeyuk (Ref. ||5^). The authors investigated the sensitivity 
of the nuclear vertex constant (which is the ANC up to a trivial kincmatical factor) for the system ^^O +n ^^O 
calculated in the microscopic GCM for ten different effective NN-potentials. All the adopted potentials, among which 
were 6 Volkov potentials, HNY, Minnesota and Gogny potentials, significantly overestimated the experimental nuclear 
vertex constant, as appears also to be the case here. For example, V2 (Volkov 2) potential overestimates the ANC 
by 50%. Hence the overestimation of the microscopically calculated ANCs has a long history and the present result 
for ^B +p ^ ^C once more seems to confirm that. 

Using the astrophysical S-factor obtained above from the experimental ANC for the energy region Ecm = — 0.8 
MeV we evaluated the reaction rate shown in Fig. 6 for the temperature range Tg = — 1, which covers the relevant 
temperature range for explosive hydrogen burning in supernovae Q. Using the expansion for the reaction rate in 
powers of Tg [Q (for the case of slowly varying astrophysical S-factors) : 

R = Na{<Jv) = Tg"^ exp ( - 4-)(^o + ^i^g^ + A2TI) (cm^s/mol) (13) 

^9 
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where Na is Avogadro's number, Tg is the temperature in 10^ K, and 

ij = 3[(!I£!|l^)2ii]i/3^4 249[M[amu]Z2z|]i/3 {lO^K)-'/^ (14) 

we found that we can reproduce the calculated temperature variation with the constants: B = 11.945 (not fitted), 
Ao — 6.64 • 10^, Ai — 8.50 • 10"* and A2 = —2.41 • lO'^. The formula fits the integrated reaction rates with an accuracy 
of 0.5% over the range Tg = — 1. 



IV. CONCLUSIONS 



We have used existing experimental cross section data for the breakup of the unstable nucleus at 285 MeV/u on 
four different targets |l5[ to determine the ANC of the radial overlap integral of the last proton orbiting its core by 
comparing them with Glauber model calculations in the eikonal approximation. Two different approaches were used 
to generate S-matrix elements. In the first approach, both G-matrix and t-niatrix effective interactions were folded 
with HF single particle densities in order to obtain the interaction potentials. In the second approach, the optical 
limit of Glauber's multiple scattering theory was used with three prescriptions for the elementary amplitudes. All 
calculations gave consistent results. There were no new parameters that were adjusted in the present calculations. 
We found the ANC for the virtual decay ^C^®B+p, and then the astrophysical factor Sis f^or the radiative proton 
capture reaction ^B(p,7)^C. This reaction gives a possible path to the hot pp chain pp-IV at high temperatures and 
away from it toward a rapid alpha process rap I at high temperatures and densities . 

Moreover, we show here a case where the indirect determination of nuclear astrophysical information from the 
breakup of unstable nuclei at intermediate energies proposed by us earlier for ^B, shows its usefulness. The method 
we use was tested so far only for the case of ®B [|l^ for which an independent knowledge of the ANC exists, and this 
is a second case we propose. We also note that the breakup cross sections are somewhat larger than those for the 
transfer reactions (around 100 mb, compared with a few mb for proton transfer) due to the fact that essentially the 
whole range of radii from about the grazing radius to infinity participates, whereas for transfer only a limited region 
around the grazing radius is involved. 

Better experimental data, including detailed momentum distributions and eventually disentangling the contribution 
of the stripping and diffraction dissociation processes would provide additional information to check the reliability of 
the models used. Also more breakup data in the same region of nuclei would provide valuable information to check 
the parameters and procedures used in these calculations and increase their overall reliability. 
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FIG. 1. The coordinate system used in the Glauber model calculations. 

FIG. 2. The stripping and diffraction dissociation contributions to the breakup probability of 285 MeV/u ^C on C, Al, Sn 
and Pb targets as a function of the proton impact parameter s. 
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FIG. 3. The asymptotic normalization coefficients determined using the eikonal approach with potentials extracted using 
JLM interactions. The results are plotted for the C, Al, Sn and Pb targets (left to right). The uncertainties are a sum in 
quadrature of the experimental errors and the uncertainties in the calculations. 

FIG. 4. Comparison of the calculated reaction cross section (full line) and total cross section (dashed line) for p+^^C (top 
panel), and of the reaction cross section of *B (middle panel) and (bottom panel) on ^^C target with experimental data. 
Data from Refs. ||l|,^^6[0,|| (see text). 

FIG. 5. Same as Fig. 3, but the results of the calculations with Franey-Love (squares), standard NN (triangles) and Ray 
(circles) interactions are added to JLM (filled circles). 



FIG. 6. The estimated reaction rate for the temperature range Tg = ~ 1. 
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TABLE I. Summary of the ANC extracted from different C breakup reactions and five of the Glauber calculations. 







JLM 


Franey-Love 


standard NN 


Ray NN 


zero-range NN 


target 


a{Sa) [mb] 


C!ff [fm-i] 




C!ff [fm-i] 


C!ff [fm-i] 




C 


48(8) 


1.132 


1.352 


1.244 


1.304 


1.446 


Al 


55(11) 


1.069 


1.167 


1.077 


1.125 


1.220 


Sn 


146(31) 


1.482 


1.448 


1.391 


1.438 


1.496 


Pb 


181(40) 


1.183 


1.140 


1.115 


1.145 


1.172 


average 




1.18 


1.26 


1.19 


1.24 


1.32 
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